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Abstract 

The non-equilibrium dynamics of a one-dimensional Ising model with uni- 
form, short-ranged three-spin interactions is investigated. It is shown that this 
model possesses an exponentially large number of metastable configurations 
that are stable against single spin flips. This glass- like situation results in a 
complete freezing of the system at low temperatures for times smaller than 
an intrinsic time-scale, which diverges exponentially with inverse tempera- 
ture. Via thermal activation the system eventually escapes from this frozen 

state, which is signals the onset of aging by domain growth. 
75.10N, 75.50L, 75.40G. 
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Glassy dynamics is one of the most fascinating subjects in modern physics [0,0. It 
manifests itself in an extremely slow relaxation exceeding laboratory time-scales caused by 
a rough energy landscape. In spin glasses this complex dynamics is due to frustration and 
randomness |j| and results in well known aging effects ||,[5] meaning a characterstic history 
dependence of dynamical observations. 

However, as is obvious from the situation in e.g. window glass, disorder is not a necessary 
ingredience for these phenomena to occur. For spin models it has been pointed out long 
time ago || that frustration without disorder might be able to produce a low temperature 
behavior reminiscent of spin glasses. Indeed, experiments with geometrically frustrated 
antiferromagnets §,[?]] showed anomalous dynamical behavior below a certain temperature 
that was interpreted as a spin glass transition. The magnetic properties of these materials 
are supposed to be adequately described by antiferromagnetically-coupled Heisenberg spins 
on a Kagome lattice and recent numerical studies of this model also yield indications for 
dynamical freezing || and possibly spin glass ordering ]10] at low temperatures. 

On the other side, also disorder without frustration can cause anomalous dynamics and 
slow relaxation: Diluted or disordered ferromagnets are prominent examples [0,0 for 
systems in which non-equilibrium dynamics is characterized by domain growth that is dras- 
tically slowed down by pinning of domain walls at vacancy-sites or weak bonds. Obviously at 
low enough temperatures it becomes very hard to discriminate such a scenario from true spin 
glass dynamics. Even non-frustrated, non-disordered models at or close to criticality will 
exhibit anomalously slow non-equilibrium relaxation including aging, originally supposed to 
be typical for spin glasses |TT| , p!3| . 

In this short note we present a very simple spin model without any frustration or disorder 
that possesses many glass like features at low temperatures. It turns out that quenching the 
system from high to low temperatures the dynamics gets stucked within microscopic time 
scales in a metastable configuration without long range order. The system remains frozen 
in this amorphous state for a macroscopic time before it will try to reach an energetically 
more favorable, ordered state via domain growth. This is very reminiscent of the structural 
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glass transition scenario ||, however this should only be taken as a pictorial analogy. 

The model which we consider is a one-dimensional Ising spin system with p-spin inter- 
actions defined by the Hamiltonian 

L 

H = — J ^2 ' Si+p-i , (Si = ±1) . (1) 

1=1 

The constant J fixes the energy scale of the spin interactions, which we set to one from now 
on, and periodic boundary conditions are imposed. The dynamics is defined to be the usual 
Glauber-dynamics fl4j] , where each spin is flipped with probability 



w(Si -> -Si) = i[l - S l tanh(/ il /T)] , (2) 

where the local field hi is defined to be one half of the energy difference between the con- 
figuration with spin Si = +1 and the same configuration with spin Si = — 1, and T is the 
temperature. We use sequential update for the analytic calculations as well as for the nu- 
merical simulations. For p = 2 one gets the well known ferromagnetic Ising chain, whose 
dynamics with random sequential update has been solved by Glauber []14|| . Note that for p 
odd the Hamiltonian (|I]) does not have the usual spin-flip symmetry Si — > —Si. Furthermore 
the local field is a sum over p terms with values +1 or —1, from which it follows that it can 
never be identical to zero in the case p odd, whereas for even p the probability for a spin to 
have zero local field is finite. Thus one has to discriminate between two different situations. 



For p odd a new situation arises, similar to p = 2 with external field [15]. For p even and 



p > 4 features of the latter case will coexist with those already known from the case p = 2 



14]. Therefore we will concentrate on the most simple case with p odd, i.e. p = 3 from now 



on and leave the study of p > 4 to a more detailed analysis fll7 ]. 



First we study the zero-temperature properties of this model. The groundstate of this 
system has a 4-fold degeneracy (in general 2 P_1 ) given by 
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(1) ••• + + + + + + + + + ••• 

< 2 > - + - + - + — - (3) 

(3) ...- + -- + -- + -... 

(4) + __ + __ + ... 

Introducing a local energy-variable T{ = Si-\SiSi + i all groundstate configurations are de- 
scribed by Ti = +1 for all sites i. Consider a configuration in r and S variables 

(r) ■•■ + + + + + + + + + ++ - + + + + + + + + + + + + + + + ••■ 

(S) ■■■ + + + + + + + + + ++ + - - + -- + -- + -- + -- + -.. (4) 

i 

which consists of two domains, both being in a minimum energy configuration, separated by 
a domain wall located at site % (note that in general one has 2 P_1 different kinds of domains, 
corresponding to the ground state degeneracy). A closer look to the transition probabilities 
(§) tells us that it costs an energy amount of 2Sjhj = 2(rj_i + Tj + Tj + i) to flip one spin, i.e. 
in configuration @ 6J for spins within the domains and 2J to flip the spins at the three 
sites i — and i + 1 surrounding the domain wall. Thus this configuration is stable against 
single spin flips and at finite, but small temperatures it needs a time 

^freeze ~ exp(2J/T) (5) 

to move the domain wall one lattice spacing to the left or right. This is the time scale (which 
is quite large for T < J) on which the particular configuration remains frozen for small 
temperatures. Moreover, it can be shown that all configurations of the type consisting 
(expressed in r- variables) of strings of arbitrary length I > 2 with r = +1 separated by 
isolated sites with r = — 1 are indeed metastable: From what has been said above it is clear 
that the criterion for metastability is r^_i + Tj + r i+ i > for all sites i. Thus in any triplet 
(ti_l, Ti, Tj+i) at most one minus sign may occur. Hence any metastable configuration can 

be represented by an arbitray sequence of two elementary units "H h" and "+". The 

total number til of all possible sequences in a system with L sites can be calculated via the 
Fibonacci-like iteration til = fiL-i + tils- The result for L — ► oo is 



, . , 3 /29 + V837 a/29- V837 1 tn ^ n 

n L ~ x L with x = \ 7 + W 7 ss 1.46557 . 6 

V 54 V 54 

Starting with a random initial state the sequential update procedure at zero temperature will 
drive the system into one of this exponentially large number of metastable states within only 
two sweeps through the whole chain. A random sequential update and small nonvanishing 
temperatures will not change this scenario significantly. Thus after 2r , where r is the 
microscopic time scale, the system will be frozen for a time rf reez e = T exp(2J/T). 

What has been said so far can be quantified by looking at the spin autocorrela- 
tion function CT(t,t w ) = l/LJ2i=i{Si(t + t w )Si(t w ))T and the time-dependent energy 
E T {t) = 1/L J2i=i( T i{~t))T where (■ ■ -)t means the expectation value with respect to the 
(time-dependent) probability distribution of spin configurations determined by the Master 
equation for the stochastic process considered here (in the limit L — > 00). These quantities 



can easily be calculated for the ferromagnetic Ising chain (i.e. p = 2) [14}|). In the present 
case such a treatment is not possible for the same reasons as in the ferromagnetic Ising 
chain in an external field or the Cayley tree with branching number two However, 
the remanent magnetization Cx(t,t w ) and the energy Ejn{t) can be calculated analytically 
for zero temperature with the tools introduced in ||15|| . In this letter we only note that an 
important ingredience for this problem to be exactly solvable is the fact that the local field 
acting on the spins never vanishes (details will be published elsewhere [|i~7f ): 

C T=0 (t,0) = 0.475 for t > 2 

C r= o(l,0) = 0.5 E T=0 (1) = -0.5 

(7) 

CV =0 (1, 1) = 0.9 E T=0 (t) = -0.6 for t > 2 

C T=0 (t,t w ) = 1 for t w > 2 
According to the above arguments the relations (0) also hold for T / as long as 
t Tfreeze- To check this we performed Monte-Carlo simulations of this model and results 
for the remanent magnetization CV(t, 0) at various temperatures are shown in figure 1. One 
observes the plateau at the value 0.475 extending to larger and larger times for decreasing 
temperatures. For low temperatures the final decay of CV(£, 0) seems to be algebraic (before 



it will cross over to an ultimately exponential decay at times comparable to the equilibration 
time, cf. |[L1]|). Fitting a straight line to this decay in a log-log plot yields an intersection 
with the (imaginary) line Ct = 0.475, by which we define the time scale t p i a teau for the 
lifetime of the metastable state. The temperature dependence of this quantity is depicted 
in the insert of figure 1 and yiels, as expected t p i a teau oc r freeze . 

We looked also for the waiting time (t w ) dependence of C{t,t w ) and found that the 
zero-temperature-predictions (0) are indeed also fulfilled for finite temperatures as long as 
t < Tf r ee Z e- Note that obviously for t < Tf reeze scaling laws like C(t,t w ) ~ c(t/t w ), which 
apply in many aging scenarios |TT| , p!8|] cannot hold. However, for t > Tf reeze this conventional 
scaling is restored. Furthermore we calculated the remanent energy Ex{t) in Monte-Carlo 
simulations. The result is depicted in figure 2. As for the remanent magnetization one 
sees the characteristic initial plateau. Furthermore by plotting the value of Ej<(t) versus 
temperature T for fixed time t one gets a characteristic nonmonotonic behavior. However, 
for t — > oo the location of the minimum approaches T = and one obtains a monotonic 
increase with temperature as expected for an equilibrium-thermodynamical internal energy. 

Inspecting again configuration (|j) and the following analysis one might expect that once 
the observation time reaches the time scale Tf reeze the domain wall will perform a random 
walk. In an arbitrary metastable configuration domain walls will randomly diffuse on a 
characteristic time scale rf reeze and annihilate when two of them meet. This scenario will 
result in a \/i growth of the domain size. To check this, we measured the average domain 
size in Monte-Carlo simulations. We define size I of a domain to be equal to the number of 
spin pairs between to succeding r- variables that have the value —1. Then, at a time t, we 
count the number riiit) of domains of size I. This defines the probability Pr(l,t) = lni(t)/L 
for a spin pair to be contained in a segment of length /. 

In figure 3 we show the result of the average domain size at time t after the quench 
d{t) = J2i IPt{1, t) in a log-log plot. Again one recognizes the frozen regime from d(t) being 
constant for t <C Tf reeze . For larger times an intermediate growth regime follows and we 
inserted a graph of d(t) oc \ft for comparision. One concludes that the above mentioned 
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picture of domain wall diffusion and annihilation is indeed to be applicable here. As soon as 
t reaches the order of the equilibration time r eq , the domains stop to grow and d(t) saturates 
at a value proportional to the equilibrium correlation length £ eq (T), which can be calculated 
analytically fl?| 



£ eq (T) = 1 1 log tanh {J/T)\- 1 oc exp(J/T) for T < J . (8) 
Apart from the prefactor 3/2 (in general p/2) the result (§) is identical to the case p = 2 Hl4 



This is a general feature of the model (|l|): although the dynamics shows drastic differences 
between p even and odd it turns out that the equilibrium behavior of static quantities is 
very similar. 

At zero temperature one again can calculate the average domain size exactly, and beyond 
that also the whole probability distribution Pr=o(M) f° r a site being contained within a 
domain of size I at time t. Here we only give the result, details of the calculation will be 



published elsewhere ]I7| 



1+2 



P T =o(M = 1) =1(1-1) (-) for l>2 (9) 



and 



Pt- 


=o(* 


= 2,t > 2) = 


Pt- 


=o(* 


= 3,t>2) = l A 


Pt- 


=o(* 


> 4,t > 2) = / (3/ 



26 \ (I 



(10) 

1+3 



Remember that after two timesteps the system is frozen at zero temperature. The average 
domain size d,T=o{t) at zero temperature is then given by: 

231 

d T=0 (t = 1) = 5 , d T=0 (t > 2) = — = 5.775 . (11) 

These analytical results are compared with data obtained from Monte-Carlo simulations in 
figure 4a. The agreement is excellent even at finite temperatures for t <^ Tf rceze . In figure 4b 
we show also results for P^(Z,t) at higher temperatures for t = 1,2, an intermediate time 
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(in the growth regime) and a time larger than the equilibration time — thus reflecting the 
equilibrium distribution. 

To conclude we have presented and analyzed a simple one-dimensional model whose non- 
equilibrium dynamics seems to share many features with a glass transition. One of them is 
for instance the complete freezing of the system in an " amorphous" state for a macroscopic 
time when cooled rapidly to low temperatures. Another is that this complex dynamics is 
achieved without putting in any disorder by hand. Of course, due to its one- dimensionality, 
it does not have a phase transition and also no particular temperature can be identified with 
a glass transition (leaving aside the question whether the latter is a true equilibrium phase 
transition or of purely dynamical origin). However, having demonstrated that even very 
simple models yield a very rich dynamical behavior, gives us some confidence that in higher 
dimensional models one might find indeed a candidate that shares more or even all features 
with a glass transition (as presently discussed in the context of geometrically frustrated 
models |TJ). 

As we have shown many new features of our model arise from the presence of mutli-spin 
interactions. Thus it seems worthwhile to have a closer look to such models in two or three 
dimensions, as already discussed in [JHJ. To support this view let us mention that it has been 
pointed out several years ago []20| that mean-field models with p-spin interactions (see also 
PHI ) show a dynamical behavior that is identical to that found in mode-coupling theories 
of the structural glass transition. Moreover, very recent work on self induced disorder in 
models with long range interactions [22]| heavily rely on multi-spin interactions, too. 

This work was performed within the SFB 341 Koln-Aachen-Julich. 
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FIGURES 

FIG. 1. The remanent magnetization for various temperatures calculated via Monte-Carlo sim- 
ulation of a system with 10 6 spins. From left to right: T = 0.14, 0.17, 0.20, 0.23, 0.26 0.29, 0.32 
and 0.35. The insert shows the temperature dependence of t p i a teau defined in the text, the straight 
line is the predicted dependency Tf re e Z e ~ exp(2J/T) (J = 1). 

FIG. 2. The energy E T (t) for various temperatures, /.From left to right: 0.17, 0.20, 0.23, 0.26 
0.29, 0.32, 0.35, 0.4, 0.5, 0.6 and 0.7. The insert shows the temperature dependence of E T (W0) 
and £ T (500). 

FIG. 3. Average domain size in dependence of the waiting time t in a log-log plot. The inter- 
mediate growth (between melting of the frozen domains and final saturation by equilibration) can 
be fitted nicely to d(t) ~ t 1 / 2 . 

FIG. 4. Probability distribution Pt(1, t) for domain sizes I at time t obtained from Monte Carlo 
simulations. Left: T = 0.17 and t = 1,2, the full lines is the analytical result. Right: T = 0.5 and 
t = 1, t = 100 and t = 30000 ~ r eq . 
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